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Abstract. The high-redshift (z>2) galaxies discovered over the last few 
years with the Lyman-break technique represent, in number density, a ma- 
jor fraction of the galaxies known in the Local Universe. Thus, understand- 
^ ■ ing the properties and the nature of these high-redshift systems is instru- 

r^ , mental to our understanding of the cosmic evolution of galaxies and their 

stellar content. I briefly review the observed characteristics of the Lyman- 
psj ■ break galaxies, relate these galaxies to their most likely low-redshift coun- 

O ! terparts, and discuss the implications of dust obscuration on the global 

properties of the Lyman-break population. Finally, the observational prop- 
erties of the high-redshift population are set in the framework of a simple 
evolutionary model for the stellar, metal and dust content of galaxies, to 
derive the intrinsic star formation history of the Universe. 
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1. Introduction 



K> ' Althought this Meeting is devoted to the discussion of the nature of the 

j^ ■ Ultraluminous Far-IR Galaxies at low and high redshift, and their role in 

galaxy evolution, I will not directly tackle this topic in my presentation. I 
will instead summarize and discuss the properties of a complementary sam- 
ple of galaxies: the high-z (z>2), UV-bright systems ([30], [31]). The number 
of known z>2 galaxies is now large enough that they can be classified as 
a population, and have been used to infer the past star formation history 
of the Universe ([22], [23]). As with all statistical studies, observational in- 
complenetesses and selection biases are a concern. In the case of the high-z 
population, volume corrections, luminosity selections and, last but poten- 
tially the most important, dust obscuration effects have been discussed by 
a number of authors (e.g., [24], [3], [27], [31], [5]). Here I will highhght the 
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impact of dust obscuration on the high-z galaxies and the inferred star for- 
mation history. I hope in this way to set a ground for comparison of the 
high-z UV-bright galaxy population with the high-z Ultraluminous Far-IR 
galaxies recently discovered with SCUBA ([28], [17], [2], [10]; see also the 
contributions of I. Small, of A. Blain, and of M. Rowan-Robinson to these 
Proceedings). 

2. The Lyman-Break Galaxies and their Low-Redshift Counter- 
parts 

Since it first was presented ( [29] ) , the Lyman-break technique had stood out 
as one of the most powerful tools for identifing high-z galaxy candidates. 
As of this writing, more than 550 candidates have been spectroscopically 
confirmed to be at z~3 over an area ~0.3 square degrees and about 50 
at z^i over an area --^0.23 square degrees ([31]). In number density, the 
bright ends of the z~3 and z~4 populations have similar values to the local 
L>L* galaxies, for a flat cosmology. Even if merging may have played a 
role in changing these values over time, just the number of stars contained 
in Lyman-break galaxies at z>2 accounts for ~20-30% of all the stars 
known today. The basic fact is that the Lyman-break galaxy populations 
are a significant fraction of the total galaxy population today. Thus, un- 
derstanding the nature of the Lyman-break galaxies remains a gateway to 
understanding the evolution of galaxies. 

The identification of high-z candidates is based on the detection of the 
Lyman break at 912 A, which is the strongest discontinuity in the stellar 
continuum of star-forming galaxies. A galaxy at, say, z=3 will have the 
Lyman break redshifted to 3648 A. If a pair of filters is chosen to straddle 
the break, the galaxy will appear extremely red in this color. In order to 
avoid as much as possible low-z interlopers, one or more filters are generally 
added longward of the Lyman break to select only candidates which are 
blue in this(these) additional color (s). With a careful selection of the color 
criteria, the Lyman-break technique is extremely successful at identifying 
high-z candidates; spectroscopic confirmations give a ~95% success rate 
for the z~3 sample and a -^80% success rate for the z~4 sample ([30], 
[8], [31]). The lower success rate at z~4 is due to the incidence of low-z 
interlopers, namely elliptical galaxies at Z'-^0.5-l whose 4,000 A break falls 
inside the selection window of the filters. 

While the determination of the intrinsic nature of the Lyman-break 
galaxies, whether they are massive systems or galaxy fragments, and what 
kind of progenitors they are, is still a source of heated debate (e.g., [21], [13], 
[1]), the identification of their observational low-z counterparts appears less 
controversial. 
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By selection, Lynian-break galaxies are UV-bright, actively star-forming 
systems, with a preferentially blue spectral energy distribution (SED). Ob- 
served star formation rates (SFRs) range from a few to 50 Mq yr~^, for 
a Salpeter Initial Mass Function (IMF) in the range O.I-IOO Mq ([8]). 
This range of values is typical of what observed in Local, UV-bright star- 
burst galaxies (e.g., [4]). The restframe UV and B-band half-light radii are 
around 0.2-0.3 arcsec, which correspond to spatial radii '--^1-3 h^Q kpc, 
depending on qo ([12], [14]). The similarity of the half-light radii at both 
UV and B suggests that the UV is a reliable tracer of the full extent of 
the light-emitting body. Ground-based optical spectra, which correspond 
to the restframe 900-1800 A range for a z~3 galaxy, show a wealth of ab- 
sorption features, and sometime P-Cygni profiles in the CIV(1550 A) line 
(cf. the figures in [30]), typical of the predominance of young, massive stars 
in the UV spectrum. Currently limited ground-based near-IR spectroscopy 
(e.g. [27]) has revealed nebular line emission in these galaxies. Hybrid line 
equivalent widths constructed using the UV fiux density f(UV) as denom- 
inator, namely, EW'(A)=F(line)/f(UV), show that the observed values for 
the high-z galaxies fall in the loci observed for local starburst galaxies ([25]). 
In summary, the observational properties of the Lyman-break galaxy popula- 
tion fully resemble, in the restframe UV-optical range, those of low-redshift, 
UV-bright starburst galaxies ([25]). 

The Lyman-break galaxies share another global characteristic with the 
Local starbursts. If we parametrize the observed UV stellar continuum with 
a power law, F(A)oc A'^, Lyman-break galaxies cover a large range of /3 val- 
ues, roughly from —3 to 0.4, namely from very blue to moderately red 
(Figure 1, left panel). This range is not very different from that covered 
by the Local, UV-bright starbursts (Figure 1, right panel). Population syn- 
thesis models (e.g. [19]) indicate that a dust-free, young starburst or con- 
stant star-formation population have invariably values of /3 < —2.0, for a 
vast range of metallicities. What does cause the UV stellar continuum of 
Lyman-break galaxies to be redder than expected for a young star-forming 
population ? 

3. Dust Reddening and Obscuration in Local Starbursts and in 
Lyman-break Galaxies 

There are two main causes for a red UV stellar continuum: 1) ageing of the 
stellar population; 2) presence of dust (variations of the intrisic IMF will 
not be discussed here). 

An ageing stellar population loses the high-mass, hot stars first, and 
then, progressively, lower-mass and colder stars. In the process, the UV 
stellar continuum becomes redder and redder and, also, the 4,000 A break 
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Figure 1. (Left Panel) The distribution of the observed UV spectral indices /3 of the z~3 
galaxies (M. Dickinson 1998, private communication). Only a small fraction of all the galaxies 
are blue enough {/3 < —2) to be classifiable as dust-free, young star-forming populations. Dif- 
ferent symbols corrsponds to different color selections. (Right Panel) The distribution of /3 for a 
UV-bright sample of Local starburst galaxies ([18]). Note that the range of /3 values is similar to 
the high-z sample. The UV spectral slopes of the local starbursts correlate with Ea(B— V), the 
color excess of the optical stellar continuum due to dust reddening ([6], [4]). Blue UV spectra 
correspond to small values of the color excess, red UV spectra to large values of Es(B— V). The 
best fit line through the data is shown. 



increases in strenght. This break spans a small wavelength range, thus is 
unaffected by dust reddening. The strenght of the 4,000 A break is there- 
fore a powerful constraint on the age of the stellar population. The Local 
starbursts can have very red UV continua (/3 >0), while still showing rather 
small 4,000 A breaks, telltales of the presence of a young stellar population 
(<10^ yr) or of constant star formation (over timescales ~10^ yr, see [4]). 
Ageing of the stellar population is not the main reason for the presence of 
a red UV SED in Local starbursts. Broad-band J, H, and K observations 
provide limited information on the strenght of the 4,000 A break in high-z 
galaxies, still accurate enough to exclude ageing as a general cause for the 
red UV spectra in this case as well (Dickinson 1997, priv. communication). 

Dust reddening is then the likely cause for red UV spectra, as demon- 
strated by the correlation between /3 and color excess (Figure 1, right panel). 
Dust reddening is generally a close-to-unsolvable problem for unresolved 
stellar populations (e.g., distant galaxies), because the effective obscura- 
tion will be a combination of dust distribution relative to the emitters, 
scattering, and environment-dependence of the extinction ([33], [6]). The 
situation gets better in the case of starbursts because the high energy en- 
vironment is generally inhospitable to dust. Shocks from supernovae can 
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destroy dust grains, while gas outflows can eject significant amounts of in- 
terstellar gas and dust from the site of star formation. If little diffuse dust 
is present within the star-forming region, the main source of opacity will 
come from the dust surrounding the region. Parametrizing the 'net' obscu- 
ration of the stellar continuum as: Fobsi^) = Fo{X) iQ-O-'iEsiB-V) fc(A)^ ^-^j^ 
Fo6s(A) and Fo(A) the observed and intrinsic fluxes, respectively, obscura- 
tion in starbursts is expressed as: 

k{X) = 2.656 (-2.310 + 1.315/A) +4.88 0.63 ^m < A < 1.60 ^m 

= 2.656 (-2.156 + 1.509/A - 0.198/A2 + 0.011/A^) + 4.88 

0.12 fim< X< 0.63 fim. (1) 

The connection between the color excess Es(B— V) and the measured spec- 
tral slope P is given by the correlation in Figure 1 (right panel). 

It is worth stressing that, although dust reddening corrections for star- 
bursts are parametrized above as a foreground dust screen, Equation 1 has 
been derived with NO assumptions on the geometrical distribution of the 
dust within the galaxies. Equation 1 is a purely empirical result ([6], [4]), 
which includes into a single expression any effect of dust geometry, scatter- 
ing, and environment-dependence of the dust composition. 

Equation 1 provides a recipe for correcting the observed SEDs for the 
effects of dust reddening. Does it fully account for the dust obscuration 
as well? In other words, does Equation 1 completely recover the light from 
the region of star formation or does it misses the flux from dust-enshrouded 
regions? The answer to these questions is a positive one: Equation 1 is able 
to recover, within a factor ~2, the UV-optical light from the entire star- 
forming region of UV-bright, i.e. moderately obscured, starbursts. 

We can prove the above statement by studying the FIR emission of the 
local starbursts. Dust emits in the Far-IR the stellar energy absorbed in the 
UV-optical. However, the Far-IR emission is not, by itself, an unambiguous 
measure of the opacity of the galaxy, as the intensity of the dust emission 
is also a function of the SFR in the galaxy. A good measure of the total 
opacity of the galaxy is instead provided by the ratio FIR/F(UV) ([25]). 
The Far-IR flux, FIR, and the UV flux, F(UV), are both proportional to 
the SFR, but their sensitivity to dust has opposite trends: roughly, FIR 
increases while F(UV) decreases for increasing amounts of dust, although 
the details of the trends are dictated by the geometrical distribution of the 
dust. In the assumption that the foreground dust screen parametrization is 
valid, the FIR/F(UV) ratio is related to the UV attenuation in magnitudes, 
A(UV), via ([25]): 

FIR/F{UV) = 1.19fl0°-^^(^^) - ll , (2) 
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Figure 2. The attenuation in magnitudes at 1,600 A, A(1600), as a function of tlie UV spectral 
index /3. Tlie data points (empty and filled circles) are from the same sample of UV-bright, local 
starbursts of Figure 1 ([18]). A(l,600) is proportional to log[F I R/ F (UV)] where FIR and F(UV) 
are the Far-IR and UV observed luminosities, respectively (see text). The filled and empty circles 
correspond to two different choices of the bolometric-to-FIR correction for the dust emission. This 
is the most uncertain parameter in the derivation of Equation 2; yet, its effect on the data is 
fairly small. The straight line across the data, marked C97, is the location of the reddening 
curve of Equation 1, with /3o = —2.1 for Es(B— V)=0 ([4]). There are no adjustable parameters 
in the positioning of the C97 line. The agreement between the data points and the predicted 
attenuations from the reddening curve (Equation 1) demonstrates that reddening corrections 
can fully recover the intrinsic UV emission from these galaxies and the amount of star formation 
buried in dust is relatively small. 



where the constant value 1.19 is the combination of the ratio of the boloniet- 
ric stellar luminosity to the UV luminosity and the ratio of the bolometric 
dust emission to the FIR emission. Since F(UV) and FIR (e.g., from IRAS) 
are measurable in galaxies, as is /3, the UV attenuation can be related to the 
UV spectral slope via Equation 2 ([25]). Figure 2 shows A(UV) measured at 
1,600 A as a function of /3 for a sample of Local starbursts. Overplot on the 
data is the trend predicted by Equation 1, with (3 related to E5(B— V) using 
Figure 1 and the limiting case /3o = —2.1 for Es(B— V)=0 ([4]). Equation 1 
and Figure 1 have no adjustable parameters. The agreement between the 
data and the predicted trend is therefore impressive, especially if we take 
into account that the latter is a recipe for reddening, and could in principle 
not account for the entire dust obscuration. Discrepancies at the low end of 
the locus of data points in Figure 2 are understandable in terms of sample 
incompletenesses . 

How does all this apply to Lyman-break galaxies? The entire purpose 
of obtaining dust obscuration corrections for the high-z galaxy sample is 
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to recover the intrinsic UV emission of the galaxies, therefore deriving a 
more meaningful UV luminosity function, a more accurate value of the 
SFR for each object (which can bear into the understanding of the na- 
ture of these objects), and, finally, the intrinsic cosmic SFR density ([22]). 
Figure 1 (right panel) shows the observed UV spectral slopes of the z~3 
galaxies. Those slopes can be 'translated' into a value of the effective color 
excess, which is calculated to have mean value E<j(B— V)~0.15 for the z~3 
galaxies, or an attenuation A(1600)~1.6 mag ([31], see also [3]). Inciden- 
tally, this mean value of E5(B— V) is similar to that observed in the local 
starburst sample ([3]); this is purely coincidental, and is borne of the fact 
that the two samples of galaxies cover similar ranges of /3. A similar mean 
value of the effective color excess has been obtained by Pettini et al. ([27]) 
from the analysis of the nebular emission lines in the NIR spectra of a 
small sample of Lyman-break galaxies. Correcting the observed UV spec- 
tra for dust attenuation increases the median SFR of ~ x5 in the z>2 
galaxies and of ~ x3 in the z<l galaxies (Figure 3). The difference in the 
correction factors at low and high-z is entirely due to the different wave- 
lengths at which the two redshift regimes are probed: -^ 1,600 A the high-z 
galaxies and ~2,800 A the lower-z galaxies. The dust correction factors 
have been 'measured' only for the high-z sample and have been assumed to 
hold unchanged for the z<l sample, modulo the change in wavelength (see 
discussion in [31]). 

4. The Evolution of the Stellar, Metal, and Dust Content of 
Galaxies 

The next question in line is whether a median attenuation of --^LG mag 
in the UV is reasonable at z> 2, when galaxies where at most a few Gyr 
old and, presumably, metal- and dust-poor. Both the Cosmic Far-IR Back- 
ground (CIB) detected by COBE ([11], [16]) and the FIR-bright galaxies 
detected by SCUBA at z>l ([28], [17], [2], [10]) demonstrate that dust 
was present at high redshift. The luminosity of the CIB is about 2.5 times 
higher than the luminosity of the UV-optical Background ([23]), implying a 
proportionally higher contribution of the redshift-integrated dust emission. 
However, neither the CIB nor the SCUBA galaxies are telling us how the 
dust content in galaxies has evolved with redshift. In the case of the SCUBA 
galaxies, the redshift and luminosity distribution and the AGN fraction of 
the sources will need to be tackled before providing such information. 

The time evolution of the UV luminosity density of galaxies and of the 
derivative SFR density (Figure 3) can be used to constrain the metal and 
dust enrichment of galaxies and, therefore, the intrinsic SFR density ([5]). 
The stars which produce the observed UV luminosity at each redshift pro- 
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Figure 3. The SFR density, in Mq yr^-"^ Mpc^^, as a function of rcdshift z. The data points are 
derived from the UV luminosity density of galaxies, converted to SFR using a Salpctcr IMF in the 
range 0.35-100 M© and continuous SF. The triangles are the observed values at 2,800 A of Lilly 
et al. ([20]) and Connolly et al. ([7]). The squares are the observed values at ^1,600 A as reported 
in Madau et al. ([23]), while the crosses are the new values from the Lyman-break galaxies by 
Steidel et al. ([31]). The circles represent the intrinsic SFR densities from the data of Lilly et al., 
Connolly et al., and Steidel et al., corrected for dust obscuration using Es(B— V)=0.15 (see [31]). 
The curves marked A and B bracket the range of solutions for the intrinsic SFR density from 
the evolution model described in Section 4. Curve B agrees well with the obscuration-corrected 
data points (circles). 



duce also metals and dust with negligible delay times, at most 100-200 IMyr 
in the case of dust ([9]). The obscuration from dust will produce an ob- 
served UV flux lower than the true flux. Once the effects of the dust on 
the observed UV emission are evaluated and removed, a new SFR density 
is calculated. The procedure is repeated iteratively till convergence ([5]). 
A number of observational contraints are used in the model: no more than 
~10% of the baryons are in galaxies; inflows/outflows keep the z=0 metal- 
licity of the gas in galaxies to about solar, with a '^15% mean residual 
gas content, and the z~2-3 metallicity to about 1/10-1/15 solar ([26]); the 
intrinsic SFR density at z=0 must be comparable with that measured from 
Ha surveys ([15]); the dust emission must reproduce the observed CIB and 
not exceed the FIR emission of local galaxies. 

These constraints are still not enough to yield a unique solution; one of 
the missing ingredients is the behavior of the SFR density at z>4, where 
there are no data points. Different assumptions will lead to different intrin- 
sic SFR histories. The range of solutions is bracketed by IVIodels A and B 
in Figure 3. Figure 4 shows, for each of the two solutions, the evolution of 
the dust column density in the average galaxy and the contribution to the 
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Figure 4. The evolution of the dust column density in galaxies (top panels), expressed as op- 
tical attenuation Ay in magnitudes, and of the contribution to the CIB at selected wavelengths 
(bottom panels) as a function of redshift, for both Model A (left panels) and Model B (right 
panels). Models A and B are the intrinsic SFR densities described in Figure 3 and Section 4. 
For both models the dust optical depth remains relatively modest at all redshifts. Yet, this is 
sufficient to fully account for the CIB luminosity. The contributing flux to the CIB is shown at 
the observer's restframe wavelengths 140 /xm, 240 /xm, 450 /^m, and 850 /jm, in arbitrary units. 
The two models predict different mean redshifts for the main contributors at each wavelength. In 
theory, this difference could be used to discriminate between the two solutions. However, some 
caution should be used as the spectral shape of the CIB is quite dependent on the dust emission 
SED adopted for the individual galaxies (in our case a single temperature blackbody combined 
with a v^ emissivity model). 

CIB at selected wavelengths as a function of redshift. The latter is however 
dependent on the assumptions about the intrisic dust emission SED, which 
is not well constrained. 

Model B resembles the SFR density derived from the obscuration cor- 
rected Lyman-break galaxies (Figure 3). This demonstrates that attenu- 
ations of about 1.6 mag in the UV are perfectly reasonable within the 
framework of a simple model of stellar and dust content evolution in galax- 
ies. 



5. The Intrinsic Star Formation History of the Universe and the 
SCUBA galaxies 

The shape of Model B resembles the SFR history expected for the 'mono- 
lithic collapse' model of galaxy evolution, although expectations for hier- 
archical galaxy formation models are not ruled out ([32], see discussion in 
[31]; see also the contribution of G. Kauffmann to these Proceedings). Thus 
the monolithic-versus-hierarchical dilemma is still unsolved by our current 
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knowledge of the SFR history. Values of the SFR density at higher redshifts 
(z>5) will be able place more definite constraints on the galaxy evolution 
scenario. 

The final question we want to ask is what fraction of the total SFR 
density the Lyman-break galaxies represent at each given redshift. And 
how much of the SF is so deeply buried in dust that its accounting is 
missing. The obscuration curve discussed in Section 3 is technically valid 
only for UV-bright star- forming galaxies; it cannot, obviously, correct for 
objects which are missing from the sample because they are too dusty. 
On the one hand. Model B is only slightly in excess of the obscuration- 
corrected SFR density calculated from the z>2 galaxies (by a negligible 
amount within the observational uncertainties), and the SFR history of 
Model B is perfectly sufficient to reproduce the observed GIB. It appears 
that the fraction of SF missed by considering the Lyman-break galaxies only 
is relatively small. On the other hand, a number of considerations invite 
to take this as a preliminary statement. We know that at low redshift a 
fraction of the star formation is deeply buried in dust, and is obscured 
even at IR wavelengths. The same could happen at high redshift, and the 
SCUBA sources seem to suggest that large dust contents are not impossible 
in high-z galaxies. The angular density of the SCUBA sources is about 1/2- 
1 of that of the z^-^S galaxies, and are spread over a (possibly) much larger 
redshift range than the Lyman-break galaxies, namely over ~5-10x larger 
cosmological volumes. The SCUBA sources are then --^5-20% of the Lyman- 
break galaxies by number density, but are forming stars with SFR~300- 
500 Mq yr-i. Thus the SCUBA galaxies could stih add ^25-100% to the 
SFR density of the obscuration- corrected Lyman-break galaxies, although 
an assessment of the AGN contribution is still missing. 

Because of their characteristics, the two populations, the UV-bright 
Lyman-break galaxies and the FIR-bright SCUBA sources, are likely to be 
complementary, rather than overlapping. At the level of current knowledge, 
it appears that about 50-80% of the SF in the early phases of the Uni- 
verse is accounted for by the obscuration-corrected Lyman-break galaxies; 
the remaining 20-50% of the SF may be contained in FIR-bright sources. 
However, more investigation of the nature, luminosity distribution and red- 
shift placement of the SCUBA sources is needed before these figures can 
be taken at face value. 
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